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PREFACE 


The  investigation  reported  herein  was  authorized  by  the  Long  Beach 
Port  Authority  under  a contract.  Agreement  No.  WES  7 6-k,  dated  12  Sep- 
tember 1975.  1 

The  investigation  was  conducted  from  September  to  March  1976  by 
personnel  of  the  Hydraulics  Laboratory,  U.  S.  Army  Engineer  Waterways 
Experiment  Station  (WES),  under  the  direction  of  Mr.  H.  B.  Simmons, 

Chief  of  the  Hydraulics  Laboratory,  Dr.  R.  W.  Whalin,  Chief  of  the 
Wave  Dynamics  Division,  Mr.  D.  D.  Davidson,  Chief  of  the  Wave  Research 
Branch,  and  Mr.  C.  E.  Chatham,  Chief  of  the  Harbor  Wave  Action  Branch. 
Mr.  J.  R.  Houston,  Research  Physicist,  conducted  the  investigation  and 
prepared  this  report.  Mr.  R.  R.  Bottin,  Jr.  aided  in  the  development 
of  the  finite  element  grids  and  the  plotting  of  the  data.  Drs.  H.  S. 
Chen  and  C.  C.  Mei  of  the  Massachusetts  Institute  of  Technology  pro- 
vided documentation  of  the  hybrid  finite  element  computer  program  they 
developed  and  materials  to  aid  in  its  utilization. 

Directors  of  WES  during  the  investigation  and  the  preparation  and 
publication  of  this  report  were  COL  G.  H.  Hilt,  CE,  and  COL  John  L. 
Cannon,  CE.  Technical  Director  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Si) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  converted 
to  metric  (Si)  units  as  follows: 


Multiply 


cubic  yards 
tons  (2000  lb,  mass) 
feet  per  second 
square  feet  per  second 


feet  per  second  per  second 


0.3048 

4046.856 

0.7645549 

907.1847 

0.3048 

0.09290304 

0.3048 


To  Obtain 


metres 

square  metres 

cubic  metres 

kilograms 

metres  per  second 

square  metres  per 
second 

metres  per  second 
per  second 
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LONG  BEACH  HARBOR  NUMERICAL  ANALYSIS  OF  HARBOR  OSCILLATIONS 


ALTERNATE  PLANS  FOR  PIER  J COMPLETION  AND  TANKER 


TERMINAL  PROJECT 


PART  I : INTRODUCTION 


1.  The  objective  of  this  study  was  to  investigate  by  use  of  a 
numerical  model  the  response  to  long-period  wave  excitation  of  three 
proposed  alternatives  for  modifications  to  Pier  J in  Long  Beach  Harbor 
for  oil  tanker  berthing  and  general  cargo  facilities.  Wave-height 
amplification  factors  and  normalized  maximum  current  velocities  were 
plotted  versus  wave  period  at  the  proposed  oil  tanker  terminals  to  as- 
certain whether  or  not  significant  harbor  oscillation  might  occur  which 
would  pose  problems  to  ship  berthing. 


2.  Two  tanker  terminal  basin  configurations  (plans)  were  con- 
sidered in  the  study;  plan  1 is  shown  in  Plate  1.  Both  300-  and  600-ft* 
openings  at  the  east  end  of  the  62- ft  deep  basin  were  tested.  Plan  2 
(Plate  2)  was  similar  to  plan  1 except  that  the  proposed  landfill  addi- 
tion to  Pier  J was  reduced  to  approximately  50  acres  with  an  open-water 
area  to  the  east  for  mooring  two  250,000-dwt  tankers.  Waves  incident 
from  a direction  normal  to  the  main  breakwater  face  with  periods  from 
30  sec  to  300-340  sec  were  considered  for  both  configurations.  Wave 
amplitude  and  current  velocities  were  calculated  every  5 sec  for  the 
period  range.  Resonant  peaks  also  were  defined  by  considering  incident 
wave  periods  in  increments  as  low  as  0.25  sec. 


A table  of  factors  for  converting  U.  S.  customary  units  of  measure' 
ment  to  metric  (Si)  units  is  presented  on  page  3. 


PART  II:  NUMERICAL  MODEL 


3.  The  response  of  the  tanker  terminal  basin  to  long  wave  excita- 
tion was  determined  by  using  a hybrid  finite  element  numerical  model  de- 
veloped recently  by  Chen  and  Mei*  at  the  Massachusetts  Institute  of  Tech- 
nology. The  model  solves  the  following  generalized  Helmholtz  equation: 

2 

V • [h(x,y)V<J>(x,y)]  + ~ <|>(x,y)  = 0 (l) 

where  <{>(x,y)**  is  the  velocity  potential  defined  by  u(x,y)  = -Vd(x,y), 
with  u(x,y)  being  a two-dimensional  velocity  vector  and  w an  angular 
frequency.  Equation  1 governs  small  amplitude  undamped  oscillation's  of 
water  in  a basin  of  arbitrary  shape  and  variable  depth  forced  by  periodic 
long  waves.  It  has  been  further  assumed  that  the  flow  is  irrotational. 

h.  The  boundary  condition  along  the  shoreline  and  along  the  de- 
tached breakwater  surrounding  the  three  proposed  berths  is  that  the 
normal  component  of  the  velocity  be  equal  to  zero.  Therefore,  the 
breakwater  is  considered  as  a solid  barrier.  No  special  boundary  con- 
ditions are  made  to  account  for  the  storage  tankers  in  plan  2 or  any 
vessels  utilizing  the  tanker  terminal  basin,  since  for  long-period 
waves  of  small  amplitude  ship  heave  equals  incident  wave  height. 
Therefore,  the  vessels  do  not  significantly  alter  the  oscillation 
characteristics  of  the  basin. 

5.  The  Helmholtz  equation: 

P v2 

V 4> (x,y ) + — <|>(x,y)  = 0 (2) 

is  the  governing  equation  for  a constant-depth  ocean  region  outside 
the  basin. 

* H.  S.  Chen  and  C.  C.  Mei,  "Oscillations  and  Wave  Forces  in  an  Off- 
shore Harbor  (Applications  of  the  Hybrid  Finite  Element  Method  to 
Water-Wave  Scattering),"  Report  No.  190,  197**,  Massachusetts  Institute 
of  Technology,  Cambridge,  Mass. 

**  For  convenience,  symbols  and  unusual  abbreviations  are  listed  and 
defined  in  the  Notation  (Appendix  A). 
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6.  For  a harbor  in  a semi -infinite  ocean  with  a straight  coast- 
line there  is  an  incident,  reflected,  and  scattered  wave.  The  scattered 
wave  has  a velocity  potential  4>s  given  by 


OO 

<f>  = y a H (kr)  cos  n0  (3) 

3 n=0  n n 

where  a are  unknown  coefficients  and  H (kr)  are  Hankel  functions 
n n 

of  the  first  kind  of  order  n . 

7.  <|>s  satisfies  the  radiation  condition  that  the  scattered  wave 

must  behave  as  an  outgoing  wave  at  infinity.  This  condition  is  known 
as  the  Sommerfeld  radiation  condition  and  may  be  expressed  mathemat- 
ically as  follows: 


- ik^  4>  = 0 (U 

r-*=°>  \9r  ) s 

8.  Chen  and  Mei  used  a calculus  of  variations  approach  and  ob- 
tained a Euler-Lagrange  formulation  of  the  boundary  value  problem.  The 
following  functional  with  the  property  that  it  is  stationary  with 
respect  to  arbitrary  first  variations  of  <J> (x ,y ) was  constructed  by 
Chen  and  Mei: 

2 

FU)  = ff  1/2  {h(V<|>)2  - — 4,2}  dA 

g 


+ 1/2  § { h(  <J>  - <j>T) 


sUd  - <t>T)} 


R YI'  3n 


da  - §{h<}> 


3 ( 4>r  - <t>T ) ) 
a 9n 


da 


94>t  3 ( 4>r  - 4>t  ) \ 

i Ch4>  t — ) da  + §{h<J)  r da 

a 9n  I 9n 


(5) 


where 

A = the  region  inside  the  harbor 
§ = line  integral 
4>r  = far  field  velocity  potential 


j 

1 
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<J>j  = velocity  potential  of  the  indicent  wave 
= unit  normal  vector  outward  from  region  A 
a = boundary  of  region  A 

<p  = total  velocity  potential  evaluated  on  boundary  a 
a. 

9-  Proof  was  given  by  Chen  and  Mei  that  the  stationarity  of  this 
functional  is  equivalent  to  the  original  boundary  value  problem. 

10.  The  integral  equation  obtained  from  extremizing  the  func- 
tional is  solved  by  utilizing  the  finite  element  method.  This  method  is 
a technique  of  numerical  approximation  that  involves  dividing  a domain 
into  a number  of  nonoverlapping  subdomains  which  are  called  elements. 

11.  The  solution  of  the  problem  is  approximated  within  each  ele- 
ment by  suitable  interpolation  functions  in  terms  of  a finite  number 

of  unknown  parameters.  These  unknown  parameters  are  the  values  of  the 
field  variable  <J>(x, y)  at  a finite  number  of  points  which  are  called 
nodes.  The  relations  for  individual  elements  are  combined  into  a 


system  of  equations  for  all  unknown  parameters. 

12.  In  the  region  outside  the  basin,  the  velocity  potentials  are 
solved  analytically  in  terms  of  unknown  coefficients.  The  region  is 
considered  a single  element  with  an  "interpolation  function"  given  by 
Equation  3.  The  infinite  series  is  terminated  at  some  finite  value  such 
that  the  addition  of  further  terms  does  not  significantly  influence  the 
calculated  values  of  <j>(x,y)  . The  resulting  equation  is  combined  with 
the  system  of  equations  for  unknown  parameters  at  nodal  points  within 
the  basin  and  this  complete  system  is  solved  using  Gaussian  elimination 
matrix  methods. 

13.  h(x,y)  is  related  to  <|>(x,y)  through  the  linearized  dynamic 
free  surface  boundary  condition 


The  horizontal  velocity  components  have  the  following  form 


15.  The  hybrid  finite  element  method  (so  named  by  Chen  and  Mei 
because  the  method  involves  the  combination  of  analytical  and  finite  ele- 
ment numerical  solutions)  is  a steady-state  solution  of  the  boundary 
value  problem.  The  response  of  a harbor  to  an  arbitrary  forcing  func- 
tion can  be  easily  determined  within  the  framework  of  a linearized 
theory. 


PART  III:  RESULTS 


16.  Plate  3 shows  the  locations  of  the  tanker  terminals  and 
Plates  4-6  show  the  finite  element  grids  for  STFP  2 (300),  STFP  2 (600), 
and  STFP  3.  Plates  7-9  are  plots  of  normalized  maximum  current  velocity 
versus  wave  period  at  each  tanker  terminal  for  the  three  modifications. 
The  current  velocity  multiplied  by  the  incident  wave  amplitude  in  feet 
gives  velocity  in  units  of  feet  per  second.  The  velocities  have  no 
vertical  component  or  variation  since  only  long  waves  are  considered  in 
this  study.  Also,  the  maximum  velocity  over  one  wave  period  is  plotted. 
Plates  10-12  are  plots  of  the  wave-height  amplification  factor  versus 
wave  period  in  seconds  at  each  tanker  terminal.  The  wave-height  ampli- 
fication factor  is  defined  at  a point  inside  the  basin  as  the  wave 
height  at  the  point  divided  by  twice  the  incident  wave  height.  This  def- 
inition of  amplification  factor  is  traditional  and  is  a result  of  the 
fact  that  the  standing  wave  height  for  a straight  coast  with  no  harbor 
would  be  twice  the  incident  wave  height  due  to  the  superposition  of  the 
incident  and  reflected  waves. 

17-  Plan  STFP  2 (300)  and  plan  STFP  2 (600)  have  normalized  maxi- 
mum current  velocity  (NMCV)  resonant  peaks  of  57  and  250  sec  and  57-5 
and  260  sec,  respectively.  The  NMCV  patterns  in  Plates  13  and  15  for 
the  peaks  of  approximately  1 min  are  very  similar  for  the  two  plans  with 
the  velocities  somewhat  larger  for  plan  STFP  2 (600).  The  NMCV  patterns 
in  Plates  lU  and  16  for  the  peaks  of  approximately  4-1/2  min  are 
virtually  identical. 

18.  Plan  STFP  3 has  NMCV  resonant  peaks  of  55,  55-5,  106,  and 
350  sec  (Plates  17-20).  The  NMCV  peak  of  approximately  1 min  is  substan- 
tially larger  than  similar  peaks  for  plans  STFP  2 (300)  and  STFP  2 (600). 
The  NMCV  patterns  in  Plates  17  and  18  are  similar  to  the  patterns  in 
Plates  13  and  l4  with  high  velocities  occurring  near  corners  of  the 
breakwater.  Large  velocities  also  occur  in  the  area  of  the  storage 
tankers.  The  106-sec  NMCV  peak  is  substantially  larger  for  plan  STFP  3 
than  the  other  two  plans.  Velocities  are  very  large  at  terminals  2 and 
3 as  a result  of  the  antinode  forming  in  the  inside  corner  of  the 
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breakwater  (see  Plate  27  of  the  contours  of  wave-height  simplification). 
The  350-sec  NMCV  peak  at  terminal  1 is  similar  in  magnitude  to  the 
l*-l/2-min  peaks  at  terminal  1 for  the  other  two  plans.  However,  the 
NMCV  peaks  at  terminals  2 and  3 are  less  than  corresponding  peaks  of 
plans  STFP  (300)  and  STFP  (600). 

19-  Wave-height  amplification  factors  are  quite  similar  for  plans 
STFP  (300)  and  STFP  (600)  Plates  10-12.  Contours  of  wave-height  ampli- 
fication presented  in  Plates  21-21*  are  also  very  similar.  Antinodes 
for  the  1-min  peak  are  located  at  the  two  inside  corners  of  the  detached 
breakwater  and  along  the  shore  near  terminal  3-  A single  antinode  in 
the  center  of  the  basin  forms  for  the  l*-l/2-min  peak. 

20.  The  1-min  peak  for  STFP  3 is  substantially  larger  than 
corresponding  peaks  for  the  other  two  plans.  Antinodes  form  in  the 
inside  corners  of  the  breakwater,  near  terminal  3,  and  along  the  corner 
of  the  landfill  Plates  25  and  26.  Contours  of  wave-height  amplification 
for  the  350-sec  peak  of  STFP  3 (Plate  28)  are  similar  to  the  contours 
for  plans  STFP  2 (300)  and  STFP  2 (600). 

21.  In  summary,  larger  1-  and  2-min  resonant  oscillations  occur 
for  plan  STFP  3 than  for  plans  STFP  2 (300)  and  STFP  2 (600).  The  U-  to 
6-min  oscillations  are  somewhat  larger  for  plans  STFP  2 (300)  and  STFP  2 
(600)  than  for  plan  STFP  3.  The  oscillation  characteristics  of  plans 
STFP  2 (300)  and  STFP  2 (600)  are  very  similar. 


PART  IV:  CONCLUSIONS 


22.  The  resonant  peaks  of  the  normalized  maximum  current  veloc- 
ity and  the  wave-height  amplification  factor  do  not  appear  to  preclude 
safe  berthing  at  the  terminal  areas,  given  proper  mooring  practices  and 
adequate  mooring  equipment.  The  resonant  response  peaks  for  incident 
waves  with  periods  of  approximately  1 to  2 min  are  significantly  larger 
for  plan  STFP  3 than  for  the  STFP  2 plans.  Whether  or  not  incident 
waves  in  this  period  range  excite  significant  ship  motion  depends  upon 
ship  characteristics  (length,  width,  and  draft)  and  mooring  conditions 
(type  of  lines,  amount  of  slack,  etc.). 

23.  The  incident  wave  spectrum  and  ship  response  versus  wave 
period  curves  must  be  known  in  addition  to  the  harbor  response  curves  of 
Plates  7-12  to  definitively  determine  whether  or  not  ship  surge  problems 
resulting  from  harbor  resonance  will  occur  in  the  tanker  terminal  basin. 

2b,  Until  detailed  information  on  moored  ship  response  as  a func- 
tion of  incident  wave  amplitude  and  frequency  is  obtained,  definitive 
conclusions  on  the  precise  amount  of  ship  motion  cannot  be  made.  How- 
ever, the  oscillations  of  the  tanker  terminal  basin  are  significantly 
smaller  than  typical  oscillations  noted  throughout  the  Los  Angeles  and 
Long  Beach  Harbors  complex  in  the  numerical  study  performed  for  the 
entire  harbor  complex.* 


* J.  R.  Houston,  "Long  Beach  Harbor  Numerical  Analysis  of  Harbor 

Oscillations;  Report  1,  Existing  Conditions  and  Proposed  Improvements" 
(in  preparation),  U.  S.  Army  Engineer  Waterways  Experiment  Station, 

CE,  Vicksburg,  Miss. 
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APPENDIX  A:  NOTATION 

Boundary  of  region  A 

Area  of  region  inside  a harbor 

2 

Acceleration  due  to  gravity,  32.2  ft/sec 
Water  depth,  ft 

Hankel  function  of  the  first  kind  of  order  n 
Imaginary  number 
Wave  number,  ft  1 
Integer 

Unit  vector  normal  to  boundary  a 
Spherical  coordinate,  ft 
Time,  sec 

Velocity  in  x-direction,  ft/sec 

Velocity  in  y-direction,  ft/sec 

Angular  velocity,  radians/sec 

Cartesian  coordinate,  ft 

Cartesian  coordinate,  ft 

Unknown  coefficient 

Gradient,  ft  1 

Wave  amplitude , ft 

Spherical  coordinate,  degrees 

2 

Total  velocity  potential,  ft  /sec 

2 

Total  velocity  potential  evaluated  on  boundary  a , ft  /sec 

2 

Incident  velocity  potential,  ft  /sec 

2 

Far  field  velocity  potential,  ft  /sec 

2 

Scattered  velocity  potential,  ft  /sec 
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